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1 Introduction 

There has been considerable progress in studying quarko¬ 
nium properties at finite temperature since the work of 
Matsui and Satz pQ. In the past the quarkonium proper¬ 
ties at finite temperature were studied in potential models 
PIPII^ (for recent works see Refs. EE) The applicability 
of the potential models at finite temperature is not ob¬ 
vious [7|. It is more appropriate to study the in-medium 
modifications of meson properties at finite temperature in 
terms of spectral functions (for recent review see mni 
HU). Meson correlators in Eucledian time have been cal¬ 
culated in Lattice QCD for a long time, but to get spectral 
functions out of them was considered impossible. However 
it was shown by Asakawa, Hatsuda and Nakahara that 
using the Maximum Entropy Method one can in princi¬ 
ple reconstruct also meson spectral functions. The method 
was successfully applied at zero temperature and 

Though systematic uncertainties in the spectral function 
calculated on lattice are not yet completely understood, it 
was shown in Ref. that a precise determination of 

the imaginary time correlator can alone provide stringent 
constraints on the spectral function at finite temperature. 


2 Lattice setup and simulation details 

We simulate both bottomonium and charmonium using 
the so called Fermilab action EU- It is a special formu¬ 
lation of the 0(a) improved Wilson action with broken 
axis-interchange symmetry: time-like and space-like coef¬ 
ficients are treated independently and are determined in a 
non-perturbative way using the dispersion relation. This 
formulation is known to significantly reduce lattice arte¬ 
facts at modest lattice spacings. As our major interest 
lies at finite temperature theory we face the problem of 
time extent being very short, which makes it difficult to 


fit the correlators. Therefore, following E1?T| we intro¬ 
duce anisotropic lattice, setting as/at = ^ yf 1. 

Such calculations require immense computational power, 
as to resolve ground state from excited state we need small 
lattice spacing. One of the ways to address this problem 
is building a dedicated Lattice QCD machine such as QC- 
DOC. This supercomputer was developed by physisists 
from Columbia University, BNL, RIKEN and UKQCD. 
Three such machines, each reaching about lOTFlops peak 
performance, are currently under construction at BNL and 
EPCC. We used for our simlations QCDOC prototypes, 
single-motherboard machines at about 50 GFlops peak. 
Such resources are still not adequate for the full QCD 
simulations, so we use quenched approximation, which 
is equivalent to neglecting quark loops. Typical statistics 
gathered was 500 to 1000 measurements, separated by 400 
updates. 

To study meson properties at finite temperature one 
considers the following correlators of some operator O 

D>{t,t') = {dit)d{t'))T, D<{t,t') = {d{t')d{t))T, (1) 

where (...)t = denotes the thermal average 

m As we cannot do any lattice simulation in Minkowski 
space we need to introduce a Wick rotation and we get 
the imaginary time correlator 

G(r) = (2) 

(where T denotes time ordering). Now we define a spectral 
function through the Fourier transform of as 

= D>iu;)-D<{u;) ^ ( 3 ) 

Ztt tt 

where Dji{uj) is the retarded correlator. With the help of 
the KMS condition on HHI we arrive at the follow¬ 

ing integral relation between the imaginary time correlator 
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and the spectral function 
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Now we can try to reconstruct the spectral function by 
calculating G(r) on the lattice; however a set of compli¬ 
cations arises: At zero temperature the kernel re¬ 

duces to simple exponential and when we consider large 
Euclidean times we only see the contribution from the low¬ 
est lying meson state in (t(w), i.e G(r) = exp(—mr). How¬ 
ever at finite temperature the time interval is limited and 
excited states are as important as the ground state. Ad¬ 
ditional problems arise in lattice calculations where cor¬ 
relators are calculated only on a discrete set of Euclidean 
times tT = k/N-r, k = 0 ,...,Nr — 1 with Nt being the 
temporal extent of the lattice. In order to reconstruct the 
spectral functions from this limited information it is neces¬ 
sary to include in the statistical analysis of the numerical 
results also prior information on the structure of cr(w) (e.g. 
such as cr(a;) > 0 for tu > 0). This can be done in many 
ways, none of them being rigorous. Therefore we use two 
of such methods, namely the Maximum Entropy Method 
(MEM) |12111 dj and constrained curve fitting These 
two methods are totally different both in their assump¬ 
tions and the procedures and we use them to cross-check 
our results. 

To do a lattice study of particles with given quantum 
numbers we must produce an appropriate set of operators 
with given symmetry properties. One such set is a local 
meson operator which is bilinear in quark-antiquark fields 


Oh{t,x) = q{T,x)rHq{T,x), 0 ^ = 1 , 75 , 7 ^, 757 ^ (5) 


for scalar, pseudoscalar, vector and axial vector channels 
correspondingly. The pseudo-scalar and vector correlators 
correspond to ground state quarkonia rjc^bi^So) and 

Si) respectively. The scalar and axial-vector chan¬ 
nels correspond to P state quarkonia, Xc,b{J = 0,1). The 
temporal correlators at finite spatial momentum p then 
take the following form 

Gh{t,p) = {OHiT,p)Ol^{T, -p)), 

0{t,p) (6) 


3 Charmonium at zero and finite temperature 

We start the discussion of the numerical results with zero 
temperature spectral functions in the vector channel. In 
Fig. El we show the spectral functions at four different 
lattice spacings = 1.9, 2.9, 4.0 and 8.2 GeV. The 
anisotropy ^ = 2 for the first three lattice spacing and 

^ Other methods of introducing prior information into the 
statistical analysis has also been discussed in Ref. m 


^ = 4 for the last one. The lattice spacing was fixed us¬ 
ing the heavy quark potential and the value of Tq = 0.5 
fm for the Sommer scale (see m for further details). As 
one can see from the figure the first peak, which corre¬ 
sponds to the J/ip state, does not move as we vary the 
lattice spacings. On the contrary, the position of the other 
peaks depends on the lattice spacing. Moreover, as we go 
to finer lattices more peaks appear. Similar results were 
obained also in other channels. Thus, all the structures 
except the first peak cannot be identified with physical 
states. It is possible that these peaks in fact belong to 
the continuum distorted by effects of finite lattice. In any 
case this problem requires further analysis which will be 
presented elsewhere m We also analyzed the spectral 
functions using constrained curve fitting by using a multi¬ 
exponential Ansatz with four or more terms. The results 
from the fits are shown in FigEl For the ground state we 
get very good agreement between MEM and constrained 
fit, while it becomes worse for higher states. Nonetheless, 
the rough agreement between the two methods gives us 
some confidence as well as indicating the size of possible 
systematic errors. 

Now we are in a position to discuss spectral functions 
at finite temperature. For = 8.2 GeV we have suffi¬ 
cient number of data points for temperatures above de¬ 
confinement to analyze the spectral function. In Fig.2 the 
spectral function in pseudo-scalar and scalar channel is 
shown. The peak corresponding to rjc state survives in 
the plasma phase; moreover, its position is essentially un¬ 
changed. This is consistent with previous findings ITHItTOl 
ED The scalar channel has been analyzed only in Ref. m 
and it was found that the Xc state dissolves in the plasma 
soon after the deconfinement temperature. The analysis of 
the scalar spectral function shown in Fig.2 confirms this 
conclusions. 


4 Bottomonium spectral functions and 
correlators 


While charmonium has been studied by other collabora¬ 
tions, the bound states of 6-quarks received much less at¬ 
tention. In fact this is the first study of the bottomonium 
at finite temperature. As in the charmonium case we start 
with zero temperature physics. The bottominium spec¬ 
trum using Fermilab action was studied in Ref. m where 
extened meson operators have been used. We have stud¬ 
ied the bottomonium spectrum using two lattice spacings 
= 8.2GeV and = 12.0GeV and ^ = 4 estimated 
from the Sommer scale^. In Fig.3 bottomonia spectral 
functions in scalar and pseudo-scalar channels are shown. 
Because of the large bottomonia mass the correlators are 
quite noisy and we found that it is considerably more dif- 
ficuilt to reconstruct the spectral functions then it was in 
the charmonia case. 

^ In Ref. 1241 the IP-IS splitting was used to determine the 
lattice spacing which leads to a largely overestimated value for 
it. 







K.Petrov: Quarkonium at zero and finite temperature 


3 




Fig. 1. The zero temperature vector spectral function for char- 
monium at different lattice spacings (top). Also shown the vec¬ 
tor spectral function at = 4GeV compared with the results 
of constrained curve fit (vertical lines). 


As the analysis of the bottomonium spectral functions 
is difhcuilt already at zero temperature at finite tempera¬ 
ture we study only the temperature dependence of the cor¬ 
relator. This method is based on the following argument 
m : in Eq2|the temperature dependence of the right hand 
side comes from two sources - the spectral function itself 
and the finite temperature kernel. Now let us introduce 
the so-called reconstructed correlator using the spectral 
function at zero temperature 


Gr 


poo 

i(t, T)= / duja{uj,T = 0) 
Jo 


cosh(a;(T — 1/(2T)) 


sinh ^ 


(7) 

If the spectral function does not depend on temperature 
then the reconstructed and directly calculated correlators 
will be equal, i.e. G{t, T)/Grecon{T, T) = 1. We show their 
ratio for the pseudo-scalar channel {S — wave) and the 
scalar channel {P — wave) of the bottomonium on Fig^ 
As one can see there is no modification of the correlator for 
the rib particle till almost twice the critical temperature. 
However, quite suprisingly, for the Xb we see a drastic 




Fig. 2. The charmonia spectral functions for ^ = 8.2 GeV at 
zero temperature and above deconfinement temperature T = 
l.lTc corresponding to Nt =24. 



Fig. 3. Bottomonium spectral functions at zero temperature, 
scalar and pseudo-scalar channels for = 12GeV. 
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change already at 1.3Tc which suggests that something 
has happened. The Xb state has approximately the same 
size and binding energy as the J/ip or ri^ state which do not 
show any change till (1.7 — 2.0)Tc. Thus we would expect 
the same for Xb- It also is possible that the large change 
in the scalar correlator is not caused by the dissolution of 
the Xb but rather by modification of the continuum part 
of the spectral function m- 


5 Conclusions and Outlook 

We analysed charmonium and bottomonium systems at 
zero and hnite temperatures using Lattice QCD with Fer- 
milab action on anisotropic lattices. We reconstructed their 
spectral functions, reliably identifying ground states for 
various values of the lattice cutoff. The S — wave states in 
both systems survive in the deconfined phase till as much 
as three times the critical temperature, while P — wave 
states may dissolve shortly after the phase transition. The 
situation with higher excited states, especially 25', is very 
questionable and requires further study. 
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